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Abstract. A brief overview is given of the Continuum Shell Model, a novel approach that extends the traditional nuclear 
shell model into the domain of unstable nuclei and nuclear reactions. While some of the theoretical aspects, such as role and 
treatment of one- and two-nucleon continuum states, are discussed more in detail, a special emphasis is made on relation to 
observed nuclear properties, including definitions of the decay widths and their relation to the cross sections, especially in the 
cases of non-exponential decay. For the chain of He isotopes we demonstrate the agreement of theoretical results with recent 
experimental data. We show how the interplay of internal collectivity and coherent coupling to continuum gives rise to the 
universal mechanism of creating pigmy giant resonances. 



CONTINUUM SHELL MODEL 

Recent progress in physics of weakly bound nuclei and other marginally stable mesoscopic systems requires an 
adequate many-body theory which would combine the bound states and continuum in a consistent description. Below 
we show how it can be done. Our formalism of continuum shell model (CSM) follows closely the projection method 
by Feshbach Q|. The many -body states of the system are separated into subspaces P and Q, and the Hamiltonian 
consists of three parts, H = Hqq + Hpp + HpQ. The internal space P contains shell model (SM) states with all particles 
on bound single -particle orbitals, for example A^-nucleon m-scheme Slater determinants, labelled as 1 1). The external 
space Q contains all states \c;E), labelled by energy £ and channel c, with some particles in the continuum. These wave 
functions asymptotically behave as outgoing spherical (Coulomb) waves with respect to one or several single-particle 
arguments. We solve the Schrtidinger equation H\a;E) =E\a;E) within the total P+ Q space where an eigenstate a 
is a superposition 

|a;£)=I>i(£)|l} + E [ dE' a c (E';E)\c;E'), (1) 

1 c J 

Elimination of external states leads to the effective energy-dependent Hamiltonian that acts within P space, 

J?(E)=H PP + A(E)-^W(E). (2) 

Here A(E) is a self-energy term coming from virtual excitations into the continuum and W(E) stands for the explicitly 
non-Hermitian term due to the open decay channels, 

(1\A(E)\2)=£P / dE' ^^W , mm2)=27t £ A\{EfAl{E). (3) 
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The amplitude of coupling to the continuum, A 1(E), is introduced as A\(E) — (c\E\Hqp\\) . For Eqs. it is assumed 
for simplicity that the P and Q spaces are orthogonal, and that for the continuum Hqq \c;E) =E\c;E), Thus, given the 
effective Hamiltonian (|2j, the internal solution can be found as 



^Jf l2 (E)a 2 (E) = ^a l (E), 
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(4) 



then, using the internal amplitudes d\ (E), the continuum part can be restored as 



a c (E';E) = Ya 1 (E) . (5) 

The energy dependence of the effective Hamiltonian © in this formalism emphasizes the non-exponential nature of 
decay processes. The experimental observables can be found through the scattering matrix that is expressed via the 
same Hamiltonian: 

S CC '(E) = «p(&) 1 5"' -£Af (E) ^Al'^jexp^), (6) 

where <i; c is a smooth scattering phase in channel c coming from features not accounted in the model Hamiltonian. 

Equation (|4} has no solutions with real eigenvalues S '. This is not surprising since stationary intrinsic states cannot 
be matched by outgoing continuum states in Eq. Q. However, solutions do exist in the complex energy plane E — > $ 
as Gamow quasistationary states 01, where resonance energies are taken as E ves = Re(<£) and widths are identified 
as r res = — 2Im((f). The related non-Hermitian eigenvalue problem is discussed in j3[0|5|. This solution for the 
resonance lies in the foundation of the alternative approach, so-called Gamow Shell Model |6], where the solution 
of the non-Hermitian eigenvalue problem finds poles of the scattering matrix @. The realization of this approach, 
although well formulated mathematically, is complicated by the difficulty of analytical continuation into the complex 
plane. In general, the direct solution of Eq. @ leads to numerous unphysical roots |7], unless cuts in the complex 
momentum space are made that separate physical and unphysical regions. The Breit-Wigner formalism gives an 
alternative definition of a resonance |8], while avoiding transition into the complex energy plain. Here both effective 
Hamiltonian and scattering matrix are evaluated only on the real energy axis. The resonance energy is determined 
from the condition that real part in the denominator of scattering matrix vanishes which is equivalent to finding E les 
as E les — Re(S'), where S — E It , s — (//2)r res is a complex eigenvalue of Jf(E ms )- There is yet another definition of a 
resonance that is closely related to the experimental cross section when, from a phase shift 8 in a given channel, the 
resonance width is defined as the phase shift derivative which at the resonance produces the maximum time delay T, 
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= 0. (7) 

E = E T rs 



Having defined the resonance region as the energy region of size Y around the resonance energy, all of the above 
definitions agree if the two conditions are fulfilled: (1) the resonance regions do not overlap; (2) the energy dependence 
of the effective Hamiltonian in each of the resonance regions is weak and can be ignored. In the realistic cases, when 
the density of states is high and these conditions are violated, one has to be specific and identify the definition of 
resonances and widths. In this work we use a Breit-Wigner definition for the width. It is also important that the CSM 
can be used for calculating the cross sections which leads to direct comparison with experiment. In that case, in contrast 
to a traditional SM, the large-scale matrix diagonalization is substituted by the matrix inversion needed for Eq. (|6j. 



APPLICATION TO WEAKLY-BOUND NUCLEI 

In the following application of the CSM we assume that the Hermitian part of the effective Hamiltonian Hpp + A 
is given by the phenomenological SM Hamiltonian and limit Q-space to one and two particles in the continuum. 
Phenomenologically adjusted interactions include corrections from virtual excitations to the continuum A. Energy 
dependence of A(E) is ignored since, unlike the non-Hermitian term W, A is not sensitive to thresholds. This choice of 
the interaction assures that, for the states below decay thresholds, coupling to continuum disappears, W = 0, and the 
results of the conventional SM are exactly reproduced. 

We assume that the mean field Hamiltonian selected here as a Woods-Saxon potential is responsible for single- 
particle widths while the additional schematic two-body interaction discussed below is introduced in Hpg to mediate 
direct two-body decay processes. With ^(e) defined as a single-particle creation operator in the continuum state with 
energy e, a one-particle channel with quantum numbers j, energy of the particle £,, and the residual nucleus in the 
state \a;N — 1), so that E = E a + £j, is asymptotically labelled as |c) = b^(Ej)\a;N — 1). The two-body channel 



can be expressed similarly, \c) = b^{e)b^i{e')\a;N — 2). The first contribution from the continuum comes through a 
single-particle decay amplitude 

Af(£ + e y ) = £4(ey)<a;tf-l|*v|l;#), (8) 

V 

where the reduced single-particle amplitudes are 

ai(e) = (0\b v (e)H PQ b]\0), (9) 

and the operators b\ create particles on the discrete SM orbitals. If the valence space is small, each single-particle state 
V is uniquely identified by spin-isospin quantum numbers j. As a result, the contribution to the effective Hamiltonian 
from one-body decay is 

(l\W(E)\2)=2nS 12 E W(£j)\ 2 \{a;N-l\bj\l;N)\ 2 , (10) 

c (open) 

where E = e ; - + E a . 

Being in general a many-body operator, W is effectively reduced to a single-particle form in some important cases. 
For example, far from thresholds one can ignore the energy dependence and use the closure to simplify the summation 
in (TO), 

(l\W(E)\2) =2n8i 2 '£\aJ\ 2 \{l;N\b}bj\l;N)\ z , (11) 

W than becomes a one-body operator that assigns a width Yj ~2n\a } \ to each unstable single-particle state j = V and 
can be combined with the shell model Hamiltonian via complex single-particle energies e v = e v — ijv/2. The single- 
particle interpretation of the width is possible when the decay process is not associated with a significant change in the 
structure from parent to daughter nucleus. The coupling amplitude a J (e), Eq. 0, is calculated numerically through a 
one-nucleon scattering problem. Within few MeV above decay thresholds, the energy-dependence of a J (e) to a good 
accuracy can be parameterized by the power law consistent with the energy scaling of the width as T(e) ~ g'+V 2 
where e is the energy above the threshold and I is the partial wave of the decay channel. 

With the one-body interaction used for single-particle decays, the admixture from the two-particle continuum can 
only appear as a second order contribution, 

uy y v;v; V E-Ep-E E-Ep-e' J ' 

that proceeds through an intermediate state j3 of the nucleus with N—l particles (sequential two-body decay). The 
contribution from this decay to matrix W is 

(l\W(E)\2)=2n £ ( dedf!8(E-E a -e-f!)A\{E)*Ai(E). (13) 
aJJ J 

The integration over single-particle energies may contain poles corresponding to open one-body decays. We con- 
centrate here on the principal value part that represents off-shell processes. Eqs. (I12> and Jl 3i allow for additional 
simplifications in the near-threshold region. The overall width behaves as Y ~ E^ +l+l where Ek is the total available 
kinetic energy. 

In contrast to sequential decay, the amplitude for the direct two-body transition is mediated by the two-body 
interaction in the HpQ part of the Hamiltonian. To describe this process, a pair amplitude is introduced, 

A\{E) = flW (e, , e 2 ) (a;N -2\p L \l;N) , (14) 

where the operator pi = {b v ^>b v i}i removes a pair coupled to quantum numbers L. Unlike in the sequential case, 
direct pair emission conserves only the quantum numbers of the pair. Near thresholds we parameterize direct decay 
amplitude a^ L > (ei , £2) with the appropriate energy dependence that reflects the phase space kinematics of the free -body 
final state. 

As an example, we consider the p-shell chain of helium isotopes from 4 He to 10 He. The internal space consists 
of two single-particle levels, /? 3 / 2 and p\/2- The SM interaction and single-particle energies are defined in l9l [Toll . 
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FIGURE 1. CSM calculation for He isotopes. For each isotope, the CSM result is on the left and experimental value is on the 
right, energy of each state or resonance is not shown but the energy scale is given on the vertical axis. For each state spin, parity 
and the decay width in MeV are displayed. 



For the one-body part of Hqq and Hpq we use the Woods-Saxon potential adjusted for single-particle states in 5 He. 
Even at several MeV above threshold, the single-particle widths for the p-wave decay from ^3/2 and pi/2 levels to a 
good accuracy can be described by the parameterizations Y3/2( e ) = 0.608e 3//2 and 71/2(6) — 0.3652e 3 / 2 (all energies 
in MeV). The sequential two-body decay is computed using Eq. (II 21 with the near-threshold approximation and the 
above assumption for the energy dependence of single-particle amplitudes. The calculation with just one-body part in 
Hpq is inadequate emphasizing the need for the direct two-body decay. This decay is introduced only for pair emission 
with total angular momentum L = 0. It is assumed that all L = neutron pairs couple to the continuum with the same 
decay amplitude a^ L=0 ^(£i , £2) parameterized as a( L=0 )(£i ,£2) = +£i\/ (3.0V27T). The numerical constant 3.0, the 
only parameter of the model, is introduced to describe the unknown strength of the residual two-body interaction 
needed for two-body decay of 6 He. With parameters of the model being identified, the full solution proceeds along 
the chain of isotopes starting from 4 He, so that at each step the relevant states in the daughter systems are known. For 
each state in an TV-particle system, the solution is iterative, energy E is used to determine possible open decay modes 
into TV — 1 and V — 2 daughters, then the non-Hermitian part W(E) is constructed, and the full effective Hamiltonian is 
diagonalized leading to a new complex energy eigenvalue. For each state in the system this procedure is repeated until 
convergence. 

Figure^shows the level scheme obtained with the CSM calculation for He isotopes along with experimental data. 
Beyo nd good overall agreement, additional fine details can be noted that are in close relation to recent experimental 
data l llllll2lfl3ll . For 6 He the 2^ state is well described by the sequential two-body decay into 4 He. However, the 
description of the broad first excited + state requires direct two-body decay, in agreement with the coherent L = 
pair-excitation structure of that state 11411 . For 7 He, the CSM confirms the peculiar nature of the 5/2~ state noticed in 
the experiment 1 1 1 ]: this state decays to the excited 2 + state rather than to the ground state of 6 He. This is not the case 
for the neighboring state l/2~ that mainly decays to the ground state. Theory and very recent experimental results 
fijll are in good agreement regarding spectroscopic factors and branching ratios. Additional CSM studies of realistic 
nuclear cases (the chain of oxygen isotopes) and related discussions can be found in Ref. 1 15]. 



The drawing schematically shows the spectrum 
of particle-hole excitations on the right. Due to 
multipole interaction, the giant collective state 
(GR) is repelled up in energy. At the same 
time, coupling to the common particle con- 
tinuum creates super-radiant (SR) collectivity 
visible through the particle scattering reaction. 
The cross section of such a reaction is sketched 
on the left. Factorized forms of the real inter- 
action and of the coupling to continuum assure 
collectivity in both GR and decay modes. The 
dynamics are determined by two multidimen- 
sional vectors d = {d„} and A = {A n }, see text. 
The multipole interaction shifts the GR accu- 
mulating the multipole strength by Q. = Ad 2 
along the real energy axis from the unperturbed 
centroid e. The continuum interaction shifts the 
SR state accumulating the decay width along 
the imaginary energy axis. 
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FIGURE 2. 



LOW ENERGY BRANCH OF GIANT RESONANCE 

The CSM presented here combines structure and reaction physics on the level that goes beyond any perturbative 
treatment; the unitarity of the scattering matrix is fulfilled automatically. The interesting effects related to the nature 
of many-body states embedded in the continuum can be noted in the above example of He isotopes. Compared to the 
conventional SM, the centroids of low-lying decaying states are shifted into continuum, see Ref. flUl : in addition, the 
coupling to the continuum affects mixing of states and their spectroscopic factors. 

It was first noted by Dicke II 1611 that strong coupling of two-level atoms to a common radiation field creates 
superradiant (SR) coherence. The phase transition to the superradiant regime was found to be a universal phenomenon 
in systems with overlapping resonances ranging in nature from molecular fl7l flil Il9l Eoll to hadronic (pentaquark 
1 2 1 ]). The universality of the SR mechanism ll22ll comes from the factorizable nature of continuum coupling in Eq. 
0. With a few continuum channels, the rank of the matrix W is much smaller than the dimension of the intrinsic 
space. Strong coupling to continuum reorganizes the system in favor of creation of few SR states while the remaining 
states become long-lived being nearly orthogonal to decay. 

Giant resonance (GR) in the response function of the nucleus to a multipole excitation appears as a result of the 
coherent coupling of particle-hole excitations. In many cases excitation energies of particle-hole configurations lie 
above the particle separation energy which may lead to another collectivization towards particle emission. We address 
this interplay in the simplified example 1231 12411 assuming the effective Hamiltonian 

^12 = £i<5i2 + ^^2-^1^2 (15) 

that contains unperturbed energies, real multipole interaction and interaction through the continuum, see Fig. Hand the 
caption for additional discussion. The interplay of the collective effects depends on the "angle" between the vectors 
d and A. In the degenerate case, e„ = £, the reaction amplitude is 

T(F) _ (£- £ -n)r+A(A-d) 2 

V ; {E-e-Cl)[E-£ + (i/2)r} + (i/2)X(A-d) 2 ' K ' 

In the limits of parallel internal and external couplings, 0=0°, the SR and GR collectivizations are coherent, and 
the experiment will reveal the shifted "Giant Dicke resonance" with full strength and full width F = A 2 . This is the case 
when particle scattering would be most efficient in excitation of the GR. The opposite case of orthogonal couplings, 
<p = 90°, will result in a dark collective state with no access to the continuum. Fig.[3]shows a numerical calculation for 
non-degenerate intrinsic states when the strength and the width are both shared between the displaced GR and the SR 
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FIGURE 3. Particle scattering cross section showing the excitation of GR for different angles between GR and SR collectivities. 



in the region of unperturbed excitation energies. With the increasing angle (j) the coupling with continuum becomes 
more fragmented and less effective in driving the GR; the missing strength develops a low-energy branch revealing the 
universal mechanism for creating the so-called pigmy giant resonance. 



FINAL NOTES AND ACKNOWLEDGMENTS 

We emphasized some of the very recent developments along the road from the nuclear structure and stable nuclei 
to nuclear reactions and continuum. We were able to highlight only a handful of examples and techniques, a tiny 
part of a huge modern-day effort. The main message is that the CSM today is no longer a theoretical concept, it is 
a powerful practical tool capable of describing many -body systems on the border of stability with an impressive quality. 
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